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Abstract 

Eta-deuteron scattering lengths are calculated. A summation of the multiple 
scattering series is carried out and the result is checked against more involved 
calculations. The necessity to go beyond the fixed nucleon approximation is 
emphasized. It is shown that a quasibound or virtual state in the 77-deuteron 
system may occur within the range of 77-nucleon scattering lengths suggested 
by other experiments. 
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I. INTRODUCTION 



Few-body interactions of the rj meson may complement our knowledge on the 77-nucleon 
interaction. Of related interest is the possibility of ?7-nuclear quasi-bound states. Such states 
have been predicted by Haider and Liu and Li et al. [Q], when it was realised that the 
^-nucleon interaction is attractive. In few-nucleon systems these states are expected to be 
narrow, and thus easier to detect. So far there has been no direct experimental verification 
of this hypothesis. On the other hand, Wilkin has suggested that an indirect effect of 
such a state is seen in the rapid slope of the pd — > r] 3 B.e amplitude detected just above the 
T) production threshold J|]. Another indication of the strong three-body rjpp correlations 
follows from the recent measurement of pp — > ppr] cross sections in the threshold region || . 

The ??-deuteron is the easiest few-body system to describe. In this letter a simple formula 
is given to provide the rjd scattering matrix at energies below the deuteron breakup. Detailed 
calculations are done for the scattering length A^, where the 77-nucleon scattering length 
a V N is considered as an input. In the limiting case of fixed nucleons, this formula is found to 
be consistent with earlier calculations of Ref. 0. However, our model includes corrections 
involving effects of the continuum in the rjpp system, which are found to be necessary. 

For large values of Re a vN in the region of 0.7 to 1.0 fm, suggested by some models, we 
find the r]d system to be close to binding. In this region becomes large and depends on 
details of the r]N interaction model. In particular one finds strong dependence on the way 
the t]N scattering matrix is extrapolated to the region below the threshold. If the actual A v d 
turns out to be sizable, it will be detected from the analysis of the final state interactions 
in the pd —>■ pdrj scattering experiment performed at Celsius recently [|7| . 

II. A FORMULA FOR THE ETA-DEUTERON SCATTERING MATRIX 

The purpose of this section is to derive a simple formula that relates the low-energy 
meson-deuteron scattering amplitude to the meson-nucleon scattering length. The former is 
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found by summation of a multiple scattering series and is expressed in terms of a few basic 
multiple scattering integrals. In order to motivate the method we recall a simple formula 
for the scattering length of a meson on a pair of fixed nucleons J8],|9| 

_ 20rM m 

^-r)d - _ a vN £, y±) 

where is the meson-nucleon scattering length and Rd is the nucleon-nucleon distance. 
Eq. (H) is obtained in a simple way by setting a boundary condition for the meson wave 
function tp at each scatterer ip'/ip = l/a V N. In the simplest version of genuinely fixed 
scatterers £ = 1, but a simple correction £ = ra^j fi V N is easy to implement. Here, the 
meson-deuteron reduced mass m^d corrects for the meson propagator (1/Rd), which has to 
be referred to the NN centre of mass system. The reduced meson-nucleon mass fi V N is 
necessary to relate the meson-nucleon scattering lengths to the meson-nucleon potentials. 

Already at this stage Eq. ([!]) is a fair representation of A v d for Re a V N of about 0.3 fm 
or less. In principle it handles also situations of A v d — > oo, i.e. the cases of meson-deuteron 
bound or virtual states close to threshold. The latter may occur already at Re a V N ~ 1 fm, 



which is close to the range of the ^-nucleon scattering lengths allowed by some models []TQ 
However, for such large Re a V N Eq. (JJ) becomes rather inaccurate. In the rest of this section 
we find necessary corrections, determine the virtual or quasibound state singularities and 
discuss other related calculations of A vd . 

Let us begin with a multiple scattering expansion that follows from the three-body 
Faddeev equations for a meson interacting with a pair of nucleons labeled 1,2. For the 
situation of meson-deuteron scattering below the deuteron breakup, the series for the rj- 
deuteron T-matrix is: 

Tr/d = t% + t% + t\G()t2 + t2Got\ + t2Got\Got2 + tlCro^CrO^l 

+ (*i + t 2 )G NN (h + 1 2 ) + (ti + t 2 )G NN (t 1 + t 2 )G NN (t 1 + 1 2 ) + ... , (2) 

where ti is a meson-nucleon scattering matrix, G is the free three-body propagator and 
Gnn = GqT^nGo is that part of the three-body propagator which contains the nucleon- 



nucleoli scattering matrix Tnn- This expansion is performed in momentum space and appro- 
priate integrations over the intermediate momenta are understood. The detailed notation 
and normalisation will be given later. Now the partial summation of the series for the 
scattering amplitude is performed. The latter is determined by an average 

< M v I T vd {E) I M v >=< T vd {E) >, (3) 

where <p d is the deuteron and iprj is the free-meson s-wave function in the relative meson- 
deuteron momentum variables. We define the 7]d scattering length as 

A vd = -(2n) 2 m vd < T vd {0) > . (4) 

A partial summation of the series (§) for < T > is obtained by 

< T° > 

< *i >= rrfe < 5 > 

where < T® d >=< t x + 1 2 > is the impulse approximation and 

ni = ( 6 ) 

1 ~ <T° > ' W 

^ 1 rid > 

This partial sum is already equivalent to formula (|1|) as it contains terms of the order 
A =< t > /R d in the denominator. The expansion parameter A does not need to be small to 
guarantee the success of Eq. (f|), which works even for |A| > 1, when the multiple scattering 
is divergent. In the rjd case |A| falls into the 0.1-0.3 range. Corrections for higher orders of A 
in the denominator of Eq. ([5]) may be obtained by comparing higher orders in Eq. (Q) with 
a series expansion of Eq. ([5]) with respect to Ei and Qi. In this way the next approximation 
is obtained 

2 < T° d > 

< Tr]d >= 1 - Sh. - Si - - (^i) 2 ] - [S 2 - - [A 2 - fiiSx] (8) 

where 



" 2 ^ • (9) 

^ to r 1 to n to 

E2 ^ ' (10) 



and 



< T° > 



A 2 = 2 izz^i^r^ — . (ii) 



This procedure may be continued into a systematic method to include higher powers of 
A in the denominator. The main advantage is a strong cancellation in the £2 — (^i) 2 term 
and also higher order S n terms, as was demonstrated in the optical model calculations of 



Ref. [|TTJ for He nuclei. As we show numerically in the next section, also for the 77-deuteron 
system there is a strong cancellation in the second (and indeed also in the higher orders) 
term of the partial sum (|J). This causes the method to converge much more rapidly than 
the direct X n series. 

Before going further we write down and discuss the basic quantities entering this formal- 
ism. Momentum variables are used everywhere. These are the momenta canonical to the 
Jacobi coordinates, q^N the relative NN momentum, p v the relative rj—NN momentum, and 
the corresponding variables for the other possible pairs like (q V N, Pn)- The normalisation is 
chosen so that a "volume" element is dpdq, the propagator Go = (E — E^N^q) — -^(p)) -1 
and the scattering matrices are t v Ni(q n Ni, q'^ni^ E — E(pn 2 ))5(pn2 — ^2)5 i- e - they conserve 
the spectator momentum. The are normalised in such a way that 

V(0,0,0) = --^- (12) 

with the standard convention Im a V N > 0. Later, a separable form t v ^ = 
v v (q ri N)a ri N(E)v v (q' nN ) is used with a Yamaguchi formfactor v v = (1 + q* N / '/cH) -1 '■ The 
NN scattering matrix Tnn(qnn, <i'nni E — E n {j> v )) is normalised with a different (standard) 
sign convention that requires 



Also here a Yamaguchi separable form is used with knn = 1-41 fm -1 and the strength fitted 
to reproduce the deuteron binding energy, the scattering length a^jy = 5.405 fm. 

Calculation of the multiple scattering integrals is straightforward although tedious and 
the formulae are lengthy. For simplicity we reproduce a few dominant quantities in the zero 
meson momentum and zero range ^-nucleon force limit, although actual calculations are 
performed also taking into account a finite force range. Then the impulse approximation 
term becomes < T® d >= — 2a rt N / ((2ti) 2 fj,^) , where 



is the scattering matrix averaged over some energy region, generated by the recoil of the 
spectator nucleon. The range of the latter is given by the Fourier transform of the deuteron 
wave function 0. In a more general nonzero meson momentum case the average is given by 
the momentum distribution of the r]N pair. 

The quantity of interest is the scattering length at threshold A^. Hence, E = —Ed, 
where Ed is the deuteron binding energy, and the energies in Eq. (|14]) extend down to the 
subthreshold region. This means an extrapolation into the unphysical region by a few MeV. 
Therefore, a model is required for this extrapolation and some possibilities are discussed 
later. In general, due to the short range of rjN forces, the absence of nearby singularities 
or an rjN quasi-bound state, the energy dependence of a^^E) in the narrow subthreshold 
region is apparently smooth. In all multiple scattering integrals of interest the average value 
a is used. In this way the dominant term Ei becomes 



where v nn is the Yamaguchi formfactor for the NN separable potential. In the low energy 
region the T/vjv matrix is dominated by the deuteron pole. For E = —Ed, the integration in 




(14) 




(15) 



with 




(16) 
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Eq. (pT5|) extends from the pole down to negative energies. When T^n is limited to the pole 
term, and the NN recoil energy E NN is neglected in Eq. ([16]), expression (|T5| ) reduces to a 
simple form 



01 



^ [ [ drdr^ d (r)—±—ti(r>) (17) 



with a clear physical interpretation. The higher order terms for E n have the same structure 
corresponding to (n + 1) scatterings on the optical potential V v d = — 2a^v<^(r)/(27r/x,jjv) at 
zero incident energy. 

Similarly to Refs. p!T|JT2]| it can be expected that the series for T® d , , , . . . converges 
so rapidly that T^ d would be precise on the 1% level even in the case of a bound state 
at threshold. Indeed, the effectiveness of this expansion is confirmed in Table I for two 
particular values of the rjN scattering lengths and an rjN formfactor allowing comparison 
with the calculation of 0. The latter one uses a rather involved set of integral equations 
for the 7] scattering on fixed nucleons corrected later for the effect of the deuteron wave 
functions, but no allowance is made for a free NNrj spectrum in the intermediate states. 
This assumption would correspond to our model with = 0. The agreement between these 
two calculations is rather good, with small differences probably due to two factors: first, 
Gaussian wave functions are used in Ref . || , while ours come naturally in the Hulthen form; 
second, the T/vjv used here contains more than the deuteron pole. 

The effect of the free three-body spectrum in the intermediate states is still missing. 
To lowest order in cl v n it is given by the Q 1 term of Eq. @. Within the average 
approximation it becomes 

- f f dg'dq M<l)M-<?) / 1h n 

01 " W J (2,n vN E NN (^) + E v (ff+<?)-E = ^ (18) 

This quantity is real below the deuteron breakup, which is the region of our interest. In a 
similar way one obtains higher order terms f2 2 etc., which are also real in this region. 
Now, the final formula to be used in the applications is presented as 

A d = 2^ ^ , (19) 
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where the second order terms are defined in analogy to Eqs. (fL5| , [18|), i.e. a 2 = 
E 2 /(4a^ iV ), uj 2 = Q 2 /a^ N and 5 2 = A 2 /(4a^ 7V ). The numerical factors 1, 2 and 4 in Eq. 
(|T9|) arise from the number of independent collisions on the successive nucleons. 

Numerical results for the scattering length are given in the next section. Before presenting 
these let us discuss the question of the unitarity of < T^E) >, when it is calculated for 
finite meson energies, but below the deuteron breakup threshold. Imaginary contributions 
to Ui arise only above the deuteron breakup, but the absorptive parts of <Tj begin already 
at E > —Ed- These are generated by the deuteron pole in Tnn(E — E v (p)) in the integral 



(jTq) for o~i and in similar formulas for the higher order <jj. At E = —Ed , the pole term 
of G T NN G contributes the ^1 term to expression (p~7|) in the coordinate 

representation. For higher energies the meson propagator becomes exp^p^r — f |)/|r — f |, 
where p v d is the incident meson momentum. Hence for small momenta we have Im o\ ~ 
ipr,dTn V d/ UnN- In addition the second term a 2 — o\ generates no terms linear in p v d as may 
be seen from Eq. ( |T?D and its second order analog. This allows us to present the low energy 
behaviour of < T vd > as 

< T vd (E) > {2n) 2 m vd = -Vj^m ~ "P^ 1 ( 2 °) 

A n d\E) 

as required by unitarity. It also permits the use of the series summation method in some 
energy region close to the threshold. 



III. RESULTS 

The formula for the meson-deuteron scattering length A n d expresses it in terms of an 
"effective 7/-nucleon scattering length" a^N, which is an average of the r]N scattering matrix 
extrapolated by a few MeV below the threshold. This is our input parameter. Inherently 
there is another parameter, the inverse range k v n in the r]N form factor, included in the 
calculated quantities multiplying powers of the a V N in the expansion ( |T§D for the inverse 
1/Ajjd- These follow from the three-body interactions, and they are given in Table II for 
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four representative values of k v n from the literature |j|, 15 jTB| , which presumably cover the 
whole range of allowed values. It follows from this table that the coefficients of the a 2 order 
are very small. The second order terms oi — o" 2 due to the intermediate NN interactions 



have been found to be small already in the optical model calculations of rjHe scattering |TT 
Higher order terms are negligible as may be seen from the results of Table I. 

As shown in Table I the static nucleon approximation is not adequate in the region of 
large Re a V N close to the critical values. The effect of the NN continuum is significant. Some 
care is needed with calculations of the corresponding scattering integrals in the three-body 
continuum U{. These converge slowly and depend strongly on the range of the rjN interaction. 
In particular, the zero range limit cannot be taken for co> 3 and higher orders. Provided an 
unphysical value n nN = 00 is not used to describe the free-spectrum contributions, the small 
weight and the cancellations in co> 2 — reduce effects of the continuum to the U\ term. A 
similar cancellation happens in the "mixed term" 62 — i0\O\ which, although small, dominates 



the second order term. All this allows the simple formula to work very effectively even 
under the demanding condition of a nearby singularity. 

In Table III one finds values of A n d for a number of rjN scattering lengths that follow from 
several analyses of the combined ttN, rjN coupled channels and from 77 photoproduction data. 
The resulting A n d may vary by an order of magnitude, reflecting a nearby rjd quasibound 
state that may arise at threshold for Re a V N close to 0.8 fm. It is also clear that the results 
depend on the extrapolation of a n ^(E) below the threshold. No detailed models are available 
in this region and two simple approaches have been attempted: 

(1) a constant a V N and 

(2) a typical low-energy form a(E) = a V N/(l — iq^N^riN) as required by unitarity. 

In the latter case the calculations are done using an imaginary q v ^ = i 0.367 fm -1 . This 

2 

value follows from the average value of the subthreshold energy argument (—Ed — 2m ^ N ) 



involved in Eq. (0), (the recoil energy amounts to some 4 MeV). Such an extrapolation 
reduces the effective values of Re a n N by 10 — 20% as compared to the threshold Re a V N 
values. The sensitivity on this extrapolation method is due to the nearby singularity. This 



singularity may represent a quasi-bound rjd state, if Re A v d < 0, or a quasi- virtual rjd state, 
if Re A vd > 0. 

A more systematic study of A v d is presented in Figs. 1 and 2, where plots of Im a v d 
and Re a v d are given for two typical values of Im a V N as functions Re a V N- These are all 
calculated for k v = 3.316 fm" 1 and for the two low energy presentations of a(E). It is 
seen that the critical value of Re a V N, when the virtual state is formed at the rjd threshold 
and crosses-over to the quasi-bound state is about 0.8 fm. That is well within the range of 
expected values. If this happens, one may observe very strong effects in the final state rjd 
interactions. These may be seen in pd — > pdrj scattering experiments, if relative rjd momenta 
as small as 10-60 MeV/c are measured. 

IV. CONCLUSION 

We have applied a multiple scattering series formalism developed earlier for the 77-helium 
system to calculate the 77-deuteron scattering length. 
The conclusions are as follows: 

a) It is possible to sum effectively the multiple scattering series for the inverse meson- 
deuteron scattering matrix below the breakup threshold. The method has the advantage 
of being able to incorporate both the influence of NN scattering and the free three-body 
intermediate continuum in a relatively simple but numerically stable way without resorting 
to exact Faddeev equations. 

b) For small values of Re a V N < 0.3 fm the fixed nucleons make a good approximation. 
For larger values of Re cl v n it is necessary to take into account the intermediate three-body 
continuum states. The sensitivity to otherwise small effects is due to the proximity of an 
rjd quasibound state. For the same reason the rjd scattering length is fairly sensitive to the 
subthreshold extrapolation of the rjN scattering matrix. 

c) A virtual or quasi-bound rjd system is likely to be formed. For Re a V N ~ 0.8 fm it occurs 
close to the threshold. Such a situation may be easily detected via final state interaction 
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studies. 
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FIGURES 

FIG. 1. The real and imaginary parts of the ry-deuteron scattering length as a function of 
Re a^Tv for the constant extrapolation of the cl^n below threshold and for k v n = 3.316 fm -1 . The 
curves are given for Im a = 0.2 fm (solid) and for Im a = 0.3 fm (dashed). 

FIG. 2. The real and imaginary parts of the ??-deuteron scattering length as a function of 
Re a V: N for the energy-dependent extrapolation cl^n^E) = 0^/(1 — ^^0^) below threshold and 
for k v n = 3.316 fm -1 . The curves are given for Im a = 0.2 fm (solid) and for Im a = 0.3 fm 
(dashed). 
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TABLES 

TABLE I. Convergence of the expansion for two values of a^N = cl^n- Units are fm and 
&r)N = 3.316fm~ . 

a vN 0.27 + i0.22 0.55 + i0.30 

A° vd 0.66 + i0.54 1.35 + i 0.74 <7 f = = 

A\ d 0.66 + i0.83 1.56 + i 1.96 o\ ^ w< = 

A 2 ^ 0.64 + i0.85 1.37 + i 2.14 <r 2 ^ w< = 

0.64 + i0.85 1.37 + i 2.14 <r 3 / u { = 

Ref. § 0.65 + i0.85 1.38 + i 2.15 

A 2 vd 0.57 + i0.97 0.61 + i 2.73 a { + ^ 



TABLE II. The calculated integral quantities in the multiple scattering sum for A^ as a 



function of k^jv. 


Units are fm 1 for (o"i, 


loi) and fm 2 for (a 2, ^2) 






k v (fm -1 ) 


00 


7.617 


3.316 


2.357 


0"1 


0.4440 


0.4331 


0.3966 


0.3679 


(72 


0.2315 


0.2245 


0.1966 


0.1588 


Wl 


0.3118 


0.2958 


0.2546 


0.2246 




0.2800 


0.1870 


0.1045 


0.0723 


5 2 


0.1759 


0.1620 


0.1249 


0.1001 


a 2 - o\ 


0.0180 


0.0075 


0.0032 


-0.0039 


UJ2 - u\ 


0.1828 


0.0874 


0.0027 


-0.0025 


5 2 — 0J\G\ 


0.0374 


0.0338 


0.0239 


0.0174 
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TABLE III. Eta-deuteron scattering lengths in fm, n n M = 3.316/m 















a(E) = a vN 


a(E) = a vN 


/(I 


— iq v Na V N) 


Bhalerao-Liu j[ 


1 




0.27 


+ 


l 0.22 


0.57 + i 0.97 


0.64 


+ i 


0.81 


"modified B-L" 


I 


15, 


l| 0.44 


+ 


i 0.30 


0.63 + i 1.93 


1.01 


+ i 


1.50 


Bennhold-Tanab 


e @ 0.25 


+ 


i 0.16 


0.66 + i 0.71 


0.66 


+ i 


0.58 


"modified B-T" 




m 


10| 0.46 


+ 


i 0.29 


0.72 + i 2.04 


1.11 


+ i 


1.54 


Abaev-Nefkens 





0.62 


+ 


i 0.30 


0.36 + i 3.36 


1.65 


+ i 


2.41 


Wilkin § 






0.30 


+ 


i .30 


0.39 + i 1.28 


0.58 


+ i 


1.11 








0.55 


+ 


i 0.30 


0.61 + i 2.73 


1.40 


+ i 


1.98 


Arima |]§ 






0.98 


+ 


i 0.37 


-2.75 + i 2.77 


-0.06 


+ 


i 6.20 


Sauerman [jl9| 






0.51 


+ 


i 0.21 


1.48 + i 2.31 


1.65 


+ i 


1.39 


Batinic 






0.888 


+ 


l 0.274 


-2.90 + i 4.12 


2.37 


+ i 


5.79 








0.876 


+ 


l 0.274 


-2.76 + i 4.24 


2.42 


+ i 


5.55 


Tiator 






0.476 


+ 


i 0.279 


0.81 + i 2.15 


1.22 


+ i 


1.56 


Krusche [13] 






0.430 


+ 


l 0.394 


0.14 + i 1.91 


0.65 


+ i 


1.73 








0.579 


+ 


l 0.399 


-0.13 + i 2.64 


0.93 


+ i 


2.41 








0.291 


+ 


l 0.360 


0.17 + i 1.35 


0.42 


+ i 


1.25 
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A(rjD) (fm) 




^ 



A(rjD) (fm) 
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